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Phenomenology from lattice QCD* 

Laurent Lellouch®' 

^■Centre de Physique Theorique Case 907, CNRS Luminy, F-13288 Marseille Cedex, France 

After a short presentation of lattice QCD and some of its current practical limitations, I review recent progress in 
applications to phenomenology. Emphasis is placed on heavy-quark masses and on hadronic weak matrix elements 
relevant for constraining the CKM unitarity triangle. The main numerical results are highlighted in boxes. 


1. Introduction 

The lattice formulation of quantum field the¬ 
ory, combined with large scale numerical simula¬ 
tions, is contributing in a variety of ways to the 
current research effort in particle physics. In the 
present talk I will focus on lattice QCD and its 
role in quantifying non-perturbative strong inter¬ 
action effects. I wish to apologize straight away 
to those colleagues whose work I will not be able 
to cover. 

The objectives of lattice QCD are numerous, 
beginning with the goal of establishing QCD as 
the theory of the strong interaction in the non- 
perturbative domain. This can be achieved by 
comparing quantities calculated on the lattice, 
such as hadronic spectra, deep inelastic scatter¬ 
ing structure functions, etc. to their experimen¬ 
tally measured counterparts. Some of these as¬ 
pects were reviewed in the parallel session talk by 
Collins [f] and I will not repeat them here. Lat¬ 
tice QCD can also be used to determine the fun¬ 
damental parameters of QCD, such as the strong 
coupling as or the quark masses. Because signifi¬ 
cant progress has been made in the evaluation of 
heavy-quark masses in the last year, I will devote 
a section to them. I refer you to the talk of Wit- 
tig at Lattice 2002 for a discussion of light-quark 
masses [2]. 

Another important objective of lattice QCD 
is to evaluate the non-perturbative, strong- 
interaction corrections to weak processes involv¬ 
ing quarks. Indeed, our inability to reliably quan¬ 
tify these long-distance effects is often the dom- 
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inant source of uncertainty in measuring fun¬ 
damental weak-interaction parameters, such as 
the elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. This problem has become all 
the more acute that the results from the CESR, 
LEP, the Tevatron and now the B factories in¬ 
dicate that deviations from the standard model, 
if present, must be small. Since this is a natural 
meeting point between lattice QCD and a signif¬ 
icant fraction of the audience, I will devote the 
largest part of my talk to the review of calcula¬ 
tions of weak matrix elements relevant for con¬ 
straining the unitarity triangle. This area is one 
where the interplay of lattice QCD and the rich 
experimental program of the next few years will 
be strong, with lattice QCD in turn contributing 
and improving thanks to the many precise mea¬ 
surements to come. In this exchange, one should 
not overlook the important role that CLEO-c will 
play, by providing the lattice community with re¬ 
sults of unprecedented accuracy, for D-meson de¬ 
cays in particular. Not only will these results 
provide stringent tests of the lattice method in 
the heavy-quark sector, but they will allow us to 
normalize our ^-physics results by the equivalent 
charmed quantities, thereby reducing many of our 
systematic errors. I will unfortunately not have 
the time to discuss D-meson processes here and 
I refer you to the reviews of [13,14] for recent 
summaries. 

The discussion up to now certainly does not 
exhaust the range of subjects studied in lattice 
QCD. One subject of great phenomenological in¬ 
terest, which is not covered here, is that of non- 
leptonic kaon decays and direct CP violation in 
these decays. While significant progress, both 
theoretical and numerical, has been made in the 
last years (see e.g. the reviews [5-7]), the re¬ 
sults have not yet reached a level of maturity 
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which permits detailed comparison with experi¬ 
ment. Another important subject is the role that 
lattice QCD plays in understanding QCD and fi¬ 
nite temperature and/or density. This aspect was 
covered by Alford in his plenary lecture at this 
conference [8]. Lattice QCD can also be used to 
make predictions for exotic hadrons, such as bbg 
(see e.g. [1]), and provides a powerful tool for un¬ 
derstanding the mechanism(s) of confinement and 
spontaneous chiral symmetry breaking. 

The remainder of the talk is organized as fol¬ 
lows. In Sec. 2, I briefly review the lattice method 
and some of its current practical limitations. In 
Sec. 3, I discuss new developments in the com¬ 
putation of the charm and bottom-quark masses. 
I then review recent progress in the calculation 
of matrix elements relevant for constraining the 
CKM unitarity triangle in Sec. 4 ^ and provide 
an outlook in Sec. 5. 

2. The challenges of lattice calculations 

Before reviewing some of the limitations of 
present day zero-temperature and zero-density 
lattice calculations, let me first say a few words 
about why lattice QCD is needed at all. As all of 
you know, quarks and gluons are confined within 
hadrons. Since confinement cannot be explained 
in perturbation theory, a non-perturbative tool 
is required to relate experiment to the underly¬ 
ing theory. Lattice gauge theory is a good candi¬ 
date, as it provides a mathematically sound defi¬ 
nition of non-perturbative QCD. In lattice QCD, 
euclidean spacetime is approximated by a discrete 
set of points. Quark fields reside on these points 
and gauge fields on the links between the points, 
as shown in Fig. 1. When the lattice is finite, this 
discretization provides both ultraviolet and in¬ 
frared cutoffs and yields a well defined path inte¬ 
gral. Because the number of degrees of freedom is 
finite and the integrand positive, this integral can 
be evaluated numerically using stochastic meth¬ 
ods. Lattice QCD thus allows the computation 
of hadronic observables directly from the QCD 
lagrangian, i.e. from first principles. The only 
source of errors are the finite statistics and the 

^Please see also the parallel session talks of D. Becirevic [3] 
and J. Simone [4]. 
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Figure 1. Two dimensional projection of a four 
dimensional cubic lattice, a is the lattice spacing. 

discretization of spacetime. The former decrease 
with the square-root of the number of samples 
generated. The latter are proportional to powers 
of aAqcD and where a is the lattice spacing 
(see Fig. I), Aqcd is the QCD scale and p^ de¬ 
notes the four-momenta of the particles studied. 
The two errors can be made arbitrarily small, by 
increasing the statistics and by reducing the lat¬ 
tice spacing. 

This is how lattice QCD works in principle. Let 
me now review some of the main difficulties that 
we encounter in practice. 

2.1. Numerical cost and quenching 

Though the number of degrees of freedom is fi¬ 
nite, it is still very large. On a 32^ x 64 lattice, 
necessary to have a box of three-and-some fermi 
with an inverse lattice spacing of 2 GeV, the num¬ 
ber of degrees of freedom is of order 10®. In addi¬ 
tion, the cost of including sea-quark effects is very 
high and increases rapidly as the sea-quark mass 
is reduced. This explains the “popularity” of the 
quenched approximation, where valence quarks 
are treated exactly and sea-quark effects are ap¬ 
proximated by a mean field. I will denote this 
approximation by Nf = 0, where Nf is the num¬ 
ber of sea quarks included in the calculation. The 
benefit of this approximation, or rather trunca¬ 
tion, is a large reduction in numerical cost since 
it saves calculating the very numerically intensive 
fermion determinant. The down side is that there 
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is no small parameter which controls the approxi¬ 
mation so that the size of corrections is not known 
a priori. Experience indicates that in many cir¬ 
cumstances the error induced does not exceed 10 
to 15%, as is shown, for instance, for the case 
of light-hadron masses [9]. It is clear, however, 
that for processes which depend dominantly on 
the existence of sea-quark effects, the quenched 
approximation will be very poor. 

While significantly more costly, more and more 
unquenched or partially quenched ^ calculations 
are being performed. Most of them are done with 
two flavors of sea quarks {Nf = 2), whose masses 
are usually substantially larger than those of the 
physical up and down quarks, instead of the three 
{u, d and s, Nf = 3) of our world. Furthermore, 
because we have much less experience with these 
calculations, we are still learning how to control 
their systematic errors. 

2.2. Light quarks and chiral extrapolations 

Present day algorithms become much less ef¬ 
fective when the mass niq of light quarks is re¬ 
duced. At the same time, unwanted hnite vol¬ 
ume effects increase. Thus at present, Nf = 2 
calculations are usually limited to niq > ms/2 
or where rus is the strange-quark 

mass. This is a regime where the “p” cannot de¬ 
cay into two “pions”. One notable exception is 
the Nf = 3 simulation by MILC where quark 
masses of about 0.2tos are reached (see [10,11], for 
instance). In any case, to reproduce the physics 
of u and d quarks, observables computed for val¬ 
ues of uiq > rriu^d must be chirally extrapolated 
to rriq = mu,d- Can this extrapolation be per¬ 
formed in a controlled fashion? The only model- 
independent theoretical guide that we have is chi¬ 
ral perturbation theory (yPT). The question then 
becomes whether the quark masses used in simu¬ 
lations are light enough for yPT at low order to 
be applicable. The question of chiral extrapola¬ 
tions will be discussed more below. It was also the 
subject of a panel discussion at Lattice 2002 [12]. 


^Partially quenched refers to calculations where valence 
and sea quarks are allowed to be different. Real world 
QCD is thus a special case of partial quenching in which 
these quarks are identified. 


2.3. Heavy quarks and lattice artefacts 

The discretization errors associated with a 
heavy quark are proportional to powers of arriQ, 
where mg is the mass of this quark. For these er¬ 
rors to remain under control, arriQ must be much 
smaller than 1, or mg <C a~^. Since present 
day inverse lattice spacings are in the range of 
2 to 4 GeV, it is clear that the b quark, with 
rrib — 5 GeV, cannot be simulated directly. A 
number of solutions to this problem exist: 

• Relativistic quarks: calculations are per¬ 
formed for a number of heavy-quark masses 
around that of the charm, where discretization 
errors are under reasonable control, and the re¬ 
sults are extrapolated in powers of 1 /toq up to 
the 6-quark mass, using heavy quark effective the¬ 
ory (HQET) as a guide. The drawback of this 
approach is that the extrapolation can be signif¬ 
icant and that discretization errors proportional 
to powers of amg may be amplihed if this extrap¬ 
olation is performed before a continuum limit is 
taken. 

• Effective theories: for heavy quarks whose 
mass mg is large compared to the other scales in 
the problem, denoted here by figcD, an expan¬ 
sion of QGD in powers of fJ.gcD/'mg can be per¬ 
formed. There are different implementations of 
this idea which go by the names of lattice HQET, 
of which the leading term is the static approx¬ 
imation, NRQGD and Fermilab approach. The 
important benefit of these approaches is that dis¬ 
cretization errors are proportional to powers of 
o-dQCD instead of powers of amg. The draw¬ 
back is that precise calculations at the physical b- 
quark mass require the calculation of corrections 
proportional to powers of l/mt, which are diffi¬ 
cult to renormalize accurately [13]. In addition, 
the number of operators whose matrix elements 
must be computed grows rapidly with the power 
of If mb. In the case of NRQGD, one is also con¬ 
fronted with the fact that the continuum limit 
cannot be taken. 

• Combination of relativistic and HQET re¬ 
sults: in this case, the extrapolations from the 
charm and infinite heavy-quark-mass regimes are 
replaced by an interpolation. This is clearly a 
very good way to reach the b quark. Discretiza- 
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tion and renormalization errors are, however, very 
different in the two theories. These interpolations 
are thus really only reliable if the results in the 
two theories are obtained in the continuum limit 
and are non-perturbatively renormalized. 

3. Heavy-quark masses 

Significant progress has been made in the de¬ 
termination of the charm and 5-quark masses in 
the last year. Unfortunately, the numerical work 
is still performed in the quenched approximation. 

3.1. The charm-quark mass 

The charm-quark mass is determined by tuning 
the bare heavy-quark mass in the calculation un¬ 
til the experimental value of an observable, such 
as the mass of the Dg meson or the spin-averaged 
mass of the 15” cc states, is reproduced. The bare 
mass is then matched onto a continuum renor¬ 
malization scheme. 

The results of two new calculations appeared 
this year. The first [15] is a modern version of the 
earlier work of [16]. While the latter made use of 
Wilson fermions, the new calculation is performed 
with O (a)-improved Wilson fermions, for which 
discretization errors are reduced from O {arric) to 
O ((arric)^). This is important here since arric ~ 
0.4. Another improvement is the use of next- 
to-next-to-next-to-leading-log (N^LL) results to 
convert the mass obtained non-perturbatively 
in the RI/MOM scheme at ~ 3GeV to its 
renormalization-group invariant (RGI) [17] and 
MS values [18,19]. This conversion had previ¬ 
ously only been performed at NLL order. The 
calculation is performed in the quenched approx¬ 
imation at a single value of ~ 2.7 GeV. 
The input used to fix the charm-quark mass is 
the Dg meson mass, with rrig determined from 
Mk and a~^ from Mk*- A number of system- 
atics are considered and the result obtained is: 

= 1.26(3)(12) GeV. 

The second calculation [20] is presumably the 
ultimate quenched nic calculation. It is also the 
first time the continuum limit has been thor¬ 
oughly studied for this quantity. The calculation 
is performed at four values of the lattice spac¬ 
ing, ranging from approximatively 2 to 4 GeV 


and makes use of O (a)-improved fermions. The 
input used is similar to the calculation of [15]. 
Five different definitions of the bare quark mass 
are considered. They are based on the vector 
Ward identity (VWI), the cs-axial Ward identity 
(AWI) and a non-singlet cc-AWI. The renormal¬ 
ized masses obtained from these definitions differ 
only by discretization errors and should there¬ 
fore agree in the continuum limit. The renor¬ 
malization is performed non-perturbatively a la 
ALPHA [21], which means that the RGI mass is 
obtained without relying on perturbation theory 
below ~ 30 GeV. The continuum extrapolation 
of the results for three of the definitions of the 
RGI mass are shown in Fig. 2, along with the 
results of other calculations. The discretization 
errors displayed by the VWI and cc-AWI results 
are of the expected size ~ (awc)^, and the RGI 
masses obtained from the three definitions ex¬ 
trapolate linearly in to surprisingly consistent 
values in the continuum limit, at a = 0. Com¬ 
parable agreement is found with the two other 
definitions. In the continuum limit the authors 
of [20] quote = 1.301 (28) (20) (7) GeV, 

using the N^LL results of [18,19] to convert the 
RGI mass to the MS scheme. The first error is 
a combination of statistical and a number of sys¬ 
tematic uncertainties; the second is due to the un¬ 
certainty in A^^^^ ]21]; the third is the difference 
MS 

of the results obtained using N'^LL and N^LL ex¬ 
pressions for the conversion from the RGI mass. 
In addition, a quenched scale ambiguity of 10% is 
found to induce a 3% uncertainty on m“®(m“®). 

Before concluding on the charm-quark mass, let 
me mention that Juge et al. [25] are performing 
a modern version of Kronfeld’s preliminary 1997 
calculation [23], where the spin-averaged mass of 
the IS' cc-states was used to tune trie and the 
lattice spacing was determined from the IP — IS 
splitting. In this new calculation, performed at 
several lattice spacings, the mass is renormalized 
perturbatively at N^LO instead of at NLO. 

My conclusion for the charm-quark mass is that 
the result of [20] is the quenched result. So I use 
it as the central value and add to it a 15% system¬ 
atic error to account for quenching uncertainties. 
Thus, I quote as a lattice summary number, 
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Figure 2. Continuum extrapolation of the RGI 
charm mass by ALPHA [20]. All quantities are 
plotted in units of the scale tq ~ 0.5 fm [22]. 
The vertical dashed line corresponds to the con¬ 
tinuum limit. The other results are: Beci’02 [15], 
Gime’98 [16] and Kron’97 [23]. Not shown is 
a result obtained from combining sumrule tech¬ 
niques with lattice results [ 24 [. 


= 1.30(4)(20)GeV . 

\_ _ _ 1 _/ 

Because of the large quenching uncertainty, this 
result does not improve on the 2002 PDG sum¬ 
mary, 1.0 GeV < < 1.4 GeV, which 

does not include lattice results [26]. It is clear 
that the error on TOc from the lattice can only be 
reduced now with an unquenched calculation. 

3.2. The 6-quark mass 

Since the 6-quark mass is the i3-meson mass up 
to a smallish correction, what we really want is 
an accurate calculation of this correction. This 
statement is all the more relevant that the 6- 
quark mass cannot be reached in large volumes on 
present day lattices. Thus we rely on the heavy- 
quark expansion of the i3-meson mass: 

Mb = mr +f + 0(%) 

= +£ + O 

where the residual mass 6m and the bare binding 
energy £ are linearly divergent individually, but 
this divergence cancels in the difference. 

What the heavy-quark expansion seeks to 
achieve is a separation of short and long-distance 
contributions to Mb- In Eq. (1), £ is meant 


to contain the long-distance components oi Mb 
and the short-distance contributions. This 

division can be made unambiguous with a hard 
cutoff, such as the lattice spacing, and £ can be 
computed in the static approximation on lattices 
which are not particularly Hne. 

Two approaches to the determination of mt 
have been pursued in the context of HQET: a 
“perturbative” one [27] and a non-perturbative 
one. The latter was proposed in the last year 
[28, 29] and represents an important theoretical 
development in the matching of HQET to QCD. 

3.2.1. The “perturbative” approach 

In this approach, the determination of mt, pro¬ 
ceeds in three steps: 

a) £ is evaluated numerically in the static ap¬ 
proximation; 

b) Sm is calculated in lattice-HQET perturba¬ 
tion theory; it is known to three loops for Nf = 0 
[30,31] and two loops for Nf 0 [32]; 

c) the 6-quark mass is then obtained through 

mrimr) = cimr){MB + 6m-£)+o(^^) , 

where c(m“®) is known to three-loop accuracy 
[33]. 

The problems with this approach are: 

• that the cancellation of power divergences be¬ 
tween £ and Sm is not complete: £ — Sm ^ odf [a 
(n=3 for N^LO calculations) which becomes large 
in the continuum limit; 

• Sm has a renormalon ambiguity of order 
Aqcd which cancels against the one in c(m“®) 
[34,35]. 

Both these problems require one to go to high 
order in perturbation theory. That this is not 
only a conceptual problem can be inferred from 
the fact that the uncertainty which the neglect of 
higher-order terms induces on is about 

200 MeV at NLO, 100 MeV at N^LO and 50 MeV 
at N^LO, for Nf = 0. 

The first problem further means that a con¬ 
tinuum limit is not possible. Thus, one is left 
with discretization errors proportional to powers 
of aAqcD and perturbative errors of order af for 
N^LO calculations. 
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3.2.2. The non-perturbative approach 

This approach was proposed very recently by 
Heitger and Sommer [28,29]. It provides a frame¬ 
work in which to match HQET and QCD fully 
non-perturbatively and in the continuum limit. In 
that sense it solves all of the problems of the per¬ 
turbative method. 

Since this development is interesting conceptu¬ 
ally, it is worth trying to get the main idea across. 
It is, however, rather technical and I will only pro¬ 
vide a sketch here: 

a) First observe that the bare quark mass 

a) = M{Lo, a) — £(Lo, a) (2) 

+0 AQCD/Tom"-, A|cd/^"“) , 

which is a parameter in the HQET lagrangian, 
knows nothing about Lq, the size of the lattice 
considered. It is only a function of the RGI 
quark mass and the lattice spacing a. In Eq. (2), 
M (Lqi a) is a quantity with mass-dimension 
one which becomes the mass M of the heavy-light 
meson containing a heavy quark of mass in 
the limit Lq —> oo and a —> 0. 

The independence of on Lq means that 

when ^ it can be studied on small 

(I/Aqcd > Lq 1/to^'^') and fine-grained 
“^1) lattices where discretization errors 
proportional to powers of ^ are small. 

b) The second step consists in equating 
^bare(^RGi^ q), (jefingd through Eq. (2), on a 
small (Lq) and large (L) box: 

M ~ £{L,a) — £{LQ,a) M{Lo, a) , (3) 

where I have assumed that L is large enough and 
a small enough for M {L, a) to be the meson 
mass M to good approximation. It is this mass 
which is going to be set to the experimental value 
of the H-meson mass to determine 

c) Zero is then inserted into the right-hand side 
of Eq. (3), in the following way: 

M ~ £{L,a)-£{Lr„a) + -£(Ti,a) (4) 

-I- £{Li,a) - £{Lo,a)M{Lo,m^’^\a) , 

^The light-quark-mass dependence of E and M is sup- 
pressed for clarity. 



Figure 3. Sequence of lattices on which the 
binding-energy differences, A£i, are computed. 

where one takes L = > Ln = 2"Lo > • • ■ > 

Lq. The trick is to note that one can consider a 
sequence of lattice spacings, a < oo < • • • < a„ <C 
1/Aqcd, and that as long as m/Li I, Eq. (4) 
is equivalent to 

M ~£{L,an)-£{Ln,an)-\ - £{Li,ai) (5) 

-h £{Li,ao) - £{Lo,ao) -k M{Lo,m'^^\a) , 

up to small discretization errors proportional to 
powers of Oi/Li and OiAgco- This breaks up the 
problem of having to accommodate long-distance 
physics on a scale L and short-distance effects as¬ 
sociated with the scale into a sequence 

of more manageable problems. Indeed, for each 
binding energy difference = £{Li+i,ai) — 
£{Li,ai), i = 0, the power divergences 

~ I/oi in £{Li+i,ai) and £{Li, oQ cancel and the 
range of scales involved, from Oi to Ti+i, is not 
too large. 

d) This means that the A£i can be computed 
at leading order in lattice HQET on the sequence 
of lattices illustrated in Fig. 3. Furthermore, in 
each step of the sequence, the continuum limit 
Oi ^ 0 can be taken. 

e) Next compute M(Lo, a) in lattice 

QCD for ~ taking a ^ 0. This 

is possible, since the range of scales involved 
a <C ^ Lq is not too large. Then, inter¬ 

polate M{Lq, 0) to the value of which 
solves 

n 

= . ( 6 ) 

i=0 

This yields -\- O (l/LgTOj'^', • • •). 

can be then converted to using N^LL 
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running [18,19] if necessary. The authors of [28] 
choose Lq ~ 0.2 fm, n = 2, L ~ 1.5fm ~ 2^Lo 
and use instead of to tune to 

The great benefit of this approach is that 
is obtained without discretization nor perturba¬ 
tive errors, in contrast to the “perturbative” ap¬ 
proach described earlier. Power corrections to the 
heavy-quark limit, however, are still present (see 
Eq. (2)). Those involving l/Lomf°^ are pecu¬ 
liar to this finite-volume approach and their size 
should be investigated carefully. 

Other approaches, based on lattice NRQCD 
and a study of upsilon [36,37] or heavy-light [38] 
spectra, have also been explored. These cal¬ 
culations were performed using NLO perturba¬ 
tion theory and are therefore less accurate than 
present day calculations. 

Having presented the various approaches used 
to determine the 6-quark mass on the lattice, I 
now review the results obtained. These results 
are summarized in Fig. 4. Only published re¬ 
sults are presented. The exception being the only 
fully non-perturbative result of [28, 29] and the 
unquenched, Nf = 2, result of [37], the latter 
because it enables an estimate of quenching un¬ 
certainties. 

Not shown are preliminary quenched, N^LO re¬ 
sults which were obtained using NRQCD and an 
extrapolation to the infinite-mass limit [1,41]. At 
this same N^LO order, the result of a prelimi¬ 
nary quenched HQET calculation was reported 
in [42]. Both these results are compatible with 
the N^LO results of [40], but have a perturbative 
uncertainty about twice as small (~ 50MeV). 

The only unquenched results available [37,40, 
41] show little change compared to the corre¬ 
sponding quenched results. However, the masses 
of the sea quarks considered are rather large. In 
[40], for instance, the calculation is performed for 
two values of the sea-quark mass, rrisea, in the 
interval TOs/ 2 to mg- These findings should be 
checked through further unquenched studies. 

Performing a simple average of the quenched 
j3 = 6.2 result of [40] and the preliminary non- 
perturbative result of [28,29] yields 

= 4.38(9)(10) GeV , 


NRQCD 

Davies 94 NLO 
Hornbostel 98 NLO 
prelim 

Ali Khan 99 NLO 
HQET 

Gimenez 97 NLO 
Marlinelli 98 NNLO 
Gimenez 00 NNLO 
Gimenez 00 NNLO 
Heitger 02 NP preli 


Figure 4. Lattice results for the b-quark mass: 
Davies 94 [36], Hornbostel 98 [37], Ali Khan 99 
[38], Gimenez 97 [39], Martinelli 98 [32], 
Gimenez 00 [40], Heitger 02 [28, 29]. Next to 
the reference is the order in perturbation the¬ 
ory used in the calculation. NP stands for non- 
perturbative. The vertical lines delimit the range 
given by the PDG [26]. 

where the first error is typical of the quenched 
calculations considered, chosen also to accommo¬ 
date the two values within one standard devia¬ 
tion, and the second is an estimate of the remain¬ 
ing quenching uncertainty, taken to be 10% of 
the Bg binding energy. The quenching error was 
reduced from the 15% used for the charm mass 
to 10% because of the evidence provided by the 
Nf = 2 results discussed above. 

4. Lattice QCD for the unitarity triangle 

Analyses of the CKM unitarity triangle (UT) 
have come to rely more and more on weak ma¬ 
trix elements calculated in lattice QCD (please 
see Stocchi’s plenary review [43] and e.g. [44,45]). 
In these analyses, the summit (p, ff) of the UT 
is constrained through Arud, the frequency of 
~ oscillations, Arud/Arug, the ratio of this 
frequency to the one for oscillations, 

ex which parametrizes indirect CP violation in 
AT ^ TTTT decays and the ratio lI4b/Ucb|, together 
with the value of sin 2(3 determined from the time- 
dependent CP asymmetry in H ^ J/ipKs. The 
theoretical expressions for Arud, Amg/Amd and 
ex involve CKM factors and short-distance Wil- 












son coefficients which are known to NLO in per¬ 
turbation theory (see e.g. [46] for a review). They 
also contain the non-perturbative QCD quantities 
Bb 

= (2° R ^ and Bk- These quantities 
are defined through the matrix elements: 

{B,\{bq)v-A{bq)v-A\B,) = ^MIjIBb, , (7) 

{^\h^i.lzq\Bq{p)) = iPfdfB, , ( 8 ) 

{K°\{sd)v-A{sd)v-A\K°} = ^MlflBK , (9) 

with q = d, s. It is for the calculation of these 
matrix elements that lattice QCD is used. 

In addition to the matrix elements mentioned 
above, lattice QCD can also be used to obtain 
the form factors relevant for B D*{D)iv and 
B TT{p)iv. Together with these form factors, 
experimental measurements of the corresponding 
differential decay rates enable independent deter¬ 
minations of \Vcb\ and |I4b|, respectively. 

4.1. Bq- Bq mixing: 

Though and Bb^ are correlated through 
the AB = 2 matrix element of Eq. (7), it is useful 
to review these quantities separately: once mea¬ 
sured, leptonic decays of i?+, which are governed 
by /s, will allow a clean determination of 
there exist many more calculations of sys- 
tematics on /b^ and Bb^ are quite different. 

In Fig. 5, I have compiled all recent results 
for Jb both quenched and unquenched. A 
similar compilation for /bs/Jb is shown in Fig. 
6. From these figures it is clear that results 
obtained with different approaches to treating 
heavy quarks (relativistic QCD vs effective theo¬ 
ries) agree fairly well, suggesting that the heavy- 
quark-mass dependence of these quantities is un¬ 
der control. 

The figures also indicate that Nf = 2 calcula¬ 
tions show an increase in Jb over its quenched 
value of up to O (15%) while these effects ap¬ 
pear to be small in /b„//b. These trends are 

^Here and below, fs = = fB^ ^nd Bb = Bb^- 

®The results of [52] include an additional systematic error 
to account for the possible effect of chiral logarithms. I 
have not included it in Figs. 5 and 6 because none of the 
other results presented consider these effects and because 
I will discuss them separately in Sec. 4.1.1. 


• N 

= 0 

A 

1 ' ' ' ' 

Nf=2 





H- 




1 

•- 

1 

1 



*' 1 1 





^+- 0 - 

H 






H 

I 1 1 

1 




• 1 ' 1 

1 

H 


N,=0 

178(20) 





150 200 250 0.9 1 1.1 1.2 1.3 


fe [MeV] 

Figure 5. Lattice results for the decay constant 
of the B meson in quenched (Nf = 0) and un¬ 
quenched QCD (Nf = 2) (left) and their ratio 
(right). Calculations are grouped by the approach 
used for the heavy quark. Only results published 
after 1996 are shown. These are: APE 00 [47], 
LL 00 [48], UKQCD 00 [49], FNAL 97 [50], CP- 
PACS 00 [51], MILC 02 [52] Ali Khan 98 [53], 
Collins 99 [54], CP-PACS 01 [55], JLQCD 02 
[65[. Also shown are my averages of the results 
for and 


confirmed by MILC who recently presented the 
first results of a calculation performed with three 
flavors of sea quarks {Nf = 3) [10, II]. Their 
preliminary hndings are fB^~^ > fB^~‘^ and 
= 1-23(4)(6) and (/b.//b)'^^=3 = 

I.18(l)tf. 

Before embarking on a discussion of system¬ 
atic errors, I wish to briefly mention that an in¬ 
teresting finite-volume technique for computing 
B-meson decay constants was presented very re¬ 
cently [56]. The first quenched results obtained 
using this method are /b = I70(II)(5)(22) MeV 
and /b, = 192(9)(5)(24)MeV [56]. 

4.1.1. fB^j'- chiral extrapolations 

In the quenched approximation, the calcula¬ 
tion of /b requires a chiral extrapolation in 
light-valence-quark mass from m^ai ^ rns/2 to 
nriu^d- In unquenched calculations, an addi¬ 
tional chiral extrapolation in sea-quark mass from 
Bisea > rBs/2 to mu,d is needed to obtain both 
/b and /b^. All of the results presented in the 
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Figure 6. As in Fig. 5, but for /bJ/b and 


figures above assume mild extrapolations. It is 
legitimate to ask whether the chiral behavior of 
Jb and Jb^ is well understood. This is a serious 
issue for UT fits, due to the importance of the 
ArUs/Amd constraint. 

According to Nf = 2 heavy meson chiral per¬ 
turbation theory (HMyPT), the light-quark-mass 
behavior of Jb, where the mass of light valence 
and sea quarks are identified {rriyai = msea), is 
given by [57] 




if 


Mss 






( 10 ) 


at leading order in the 1/Mb expansion and up to 
NNLO chiral corrections. In Eq. (10), is the 
value of Jb^/AIb in the combined heavy-quark 
and chiral limits, g is the B*Bn coupling in this 
same limit. Mss is the mass of a pseudoscalar me¬ 
son composed of a quark and anti-quark of mass 
iTisea, f is the pion decay constant in the chi¬ 
ral limit and is a scale which parametrizes 
the unknown coupling of the contact term pro¬ 
portional to Mgg. ^ For /b^ and Nf = 2, one 
has to resort to partially quenched HMyPT, since 
the valence and sea light quarks are not degener¬ 
ate. The prediction for the sea-quark-mass de¬ 
pendence of is [58] {rriyai fixed to rris) 


AO) 


= l-(l + 3/) 


Mvs 

Anf 


In 


Mt 


Aj 

/b. 


( 11 ) 


^At the order in which Eq. (10) is obtained, g cannot be 
distinguished from the physical B* Bit coupling, gi,, and / 
from /t, or even fss- 


at leading order in the 1/Mb expansion and up to 
NNLO chiral corrections. In Eq. (11), Mys is the 
mass of a pseudoscalar meson composed of quarks 
of mass rriyai = rris and rrisea- The dependence 
on sea-quark mass is recovered through the LO 
yPT relation, Myg = (M|g -|- Myy)/2. This 
dependence is much milder for /b^ than it is for 
/b, because Mys is much less sensitive to rrisea 
than is Mss- 

CLEO has recently determined gc = 0.59(7) 
using its measurement of the D*'^ width, gc is the 
D* Dtt coupling and differs from g by chiral and 
1/Mb corrections. This value is consistent with a 
recent lattice calculation which further finds little 
dependence of the D* Dn coupling on the charm 
quark mass [60]. However, agreement amongst 
different determinations of g is generally rather 
poor, the spread of which lies in the interval 0.2 < 
g < 0.7 [61-63]. 

JLQCD have achieved an impressive statisti¬ 
cal accuracy for their Nf = 2 determinations 
of /b and /b^ [64]. This has enabled them to 
study rather extensively the dependence of these 
quantities on light-quark mass. ® These prelim¬ 
inary studies are illustrated in Fig. 7, where the 
quantities Jb^/Mb and /b, \/Mb^ are plotted 
against Mgg. In the quenched case, where no 
sea quarks are present, Jb^^i '/Mb////// is a func¬ 
tion of Myy. In fact, the near linearity of these 
quenched results, observed in all such calcula¬ 
tions, is one of the reasons why chiral extrapola¬ 
tions have been done assuming mild light-quark- 
mass dependence. 

JLQCD perform a number of fits to their re¬ 
sults for JbVMb and /b, \/Mb^ versus Mgg, in¬ 
cluding fits to the HQyPT expressions of Eqs. 
(10) and (11) [64]. ® For /b, ^Mb,, Fig. 7 indi- 


®A direct comparison of the JLQCD results to those of 
MILC [52], which have larger statistical errors, is difficult: 
the two sets of results are obtained at fixed bare coupling, 
with different quark actions, in different volumes, with 
different procedures for setting the scale, in addition to 
which the masses of the b and s quarks in the JLQCD 
calculation are not yet precisely tuned to their physical 
values. However, a rough comparison indicates that the 
JLQCD and MILC results display similar features in the 
region of light-quark mass where they overlap. 

®What is actually plotted in Fig. 7, for Nf=2^ is 
'fjs [roMss)'^-, with the replacement 
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- N,=0 and 2 [quadratic] 


N,=2 [log (f=131MeV. g=0.59)]; N,=0 [linear] 



dence of on light- 

quark (valence-quark) mass in Nf = 2 (Nf = 0) 
NRQCD [64] ( [69]). M = Mss for Nf = 2 and 
M = Mvv when Nf = 0. The curves denoted log 
correspond to fits to Eqs. (10) and (11). 

cates that the extrapolated value at = Mf 
is affected only mildly by the presence of chi¬ 
ral logs. The situation is radically different for 
fB\/MB where the extrapolated value depends 
critically on the coupling g, becoming lower for 
larger values of g. For g = 0.6, JLQCD find 
that the extrapolated value of fs^/Ms is 17% 
lower than that obtained through a quadratic fit 
in M'gg [65]. They conclude that the inclusion 
of logarithms in these chiral extrapolations could 
significantly lower the value of fs and increase 
that of fBs/fB- Similar conclusions were reached 
about fsJfB by the authors of [66], using the 
predicted iVy = 3 chiral logs and information 
from lattice calculations. Based on these observa¬ 
tions Yamada, who reviewed heavy-quark physics 
at Lattice 2002, chose to quote for fs a central 
value taken from linear/quadratic chiral extrap¬ 
olations, adding to it a —17% systematic error 

ro —> 2.47 GeV~^, as obtained by rescaling ro by a long¬ 
distance observable close to the chiral limit. The physical 
units are therefore only meant to be indicative. The as¬ 
sumption behind fitting i thus obtained, to 

the HQyPT expressions of Eqs. (10) and (11) is that ro 
does not have a sea-quark-mass dependence which could 
significantly modify the coefficient of the logs in these 
equations. This is not unreasonable since the scale ro is 
defined through the force between static quarks at ro and 
is found to be ~ 0.5 fm [22]. At such distances, this force 
should not be very sensitive to pions. 


to account for possible effects of chiral logs [65]. 
Given the size of this error and its effect on UT 
fits, I wish to make a number of comments on the 
JLQCD analyses: 

• The chiral log extrapolation of the results of 
JLQCD for fB\/M b is problematic from a numer¬ 
ical analysis point of view in the sense that this 
functional form varies most rapidly in a region 
unconstrained by data. 

• The size of the NLO chiral correction of Eq. 
(10) is O (60%) at the heaviest point. With such 
a sizeable correction, higher-order terms cannot 
be ignored. Moreover, xPT is not expected to 
hold up to 1 GeV since the expansion is in Mgg 
over Stt'^P ~ (1.1 GeV)^. 

• Another important piece of information is the 
behavior of fj^ vs Mf obtained by JLQCD with 
the same gauge configurations [64]. While the 
coefficient of the chiral log in the xPT expression 
for /,r [67] may even be larger than the one for 
/bVA/b (2 vs the 1.6 obtained for g = 0.6), their 
nearly linear results for f^^ are inconsistent with 
NLO chiral behavior {x^/dof = 5.1) [64] in the 
range of sea quark masses where they claim that 
logarithmic behavior may have been observed for 
/bVMb. 

• Because JLQCD work at fixed bare gauge 
coupling and with a fixed number of lattice sites, 
the spacing and the volume of their lattice shrink 
by roughly 25% as Mgg varies from ~ 1 GeV^ to 
~ 0.3 GeV^. This rather large variation will am¬ 
plify the unphysical dependence on induced 
by mass-dependent finite-volume and discretiza¬ 
tion errors. 

• If instead of requiring that the NLO behavior 
of Eq. (10) be applicable all the way up to Mgg ~ 
1 GeV^, one models the possible effect of higher- 
order chiral terms, allowing or not the logarithmic 
behavior to set in at lower masses, one finds that 
the downward shift in JLQCD’s results for /s due 
to chiral logs is generically less than about 10% 
(for g < 0.6). For a similar analysis gives 
variations which are negligible compared to other 
systematics. 

It is important to note, also, that the uncer¬ 
tainty due to chiral logs depends on the details 
of the calculation and, in particular, on whether 
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the lattice spacing was set with a quantity which 
is sensitive to chiral logs and whether these loga¬ 
rithms were taken into account. This means that 
one cannot generically ascribe a chiral-log error 
to quantities such as Jb- However, the quenched 
results reviewed above, which form the basis of 
the summary numbers presented below, include 
an error associated with scale setting and agree 
within errors despite the fact that different groups 
set the scale with quantities which have varying 
degrees of sensitivity to chiral logs. This indi¬ 
cates that the chiral-log uncertainty due to scale 
setting is to some extent already accounted for. 
The same appears to be true for the Nf = 2 over 
Nf = 0 ratios which are also used to obtain the 
final summary numbers. Thus, 1 will assume that 
the discussion surrounding the chiral behavior of 
the JLQCD results can be used to obtain a rough 
estimate of the range of deviations that could be 
induced by the chiral logs intrinsic to fs and Jb^ ■ 
Given that discussion, the exploratory nature of 
the investigations of chiral-log effects, it seems 
reasonable to fix the central values of Jb and /b„ 
from linear/quadratic chiral extrapolations and 
to add a —10% systematic error to Jb to account 
for the uncertainty in the chiral extrapolation and 
none to /b„ . 

To reduce the chiral-log error in future analy¬ 
ses, the obvious solution is to extend the lattice 
results to smaller masses to match unambiguously 
onto NLO chiral behavior. This is numerically 
very costly, however. In the meantime, it would 
be helpful to chirally extrapolate dimensionless 
ratios of quantities in which the chiral logs can¬ 
cel or nearly cancel. For such ratios, the uncer¬ 
tainties due to contributions of chiral logarithms 
at small light-quark mass will be much reduced. 
For /b in Nf=2 calculations, and g around 0.6, 
a good candidate would be /s //,r (see discussion 
above and also [52]). To make such analyses 
reliable, it is important to treat chiral logarithms 
consistently throughout the calculation [52], not 
to consider light quarks which are too massive, 
to reduce the uncertainty on g and to determine 

'^'^Nonetheless, MILC are making impressive progress in 
that direction [11]. 

"At leading order in the IjMg expansion, fg and 
JbVM g have identical chiral expansions. 


mg-rni, 



mi [GeV"] 

Figure 8. Dependence of Bb and Bb, on 
light-quark (valence-quark) mass in Nf = 2 

(Nf = 0) NRQCD [69]. The light-quark mass 
is parametrized by the mass squared of the pseu¬ 
doscalar meson composed of quarks with that 
quark mass. The different curves correspond to 
linear (dotted) and quadratic (solid) fits. 

the extent to which 1/Mb corrections can modify 
the prediction of Eq. (10). It is also important to 
verify that variations in lattice spacing and vol¬ 
ume such as those mentioned above do not distort 
dependences on light-quark mass significantly. 

4.2. Bq - Bq-mixing: Bb^ 

There are many fewer calculations of Bb^ than 
there are of . However, the situation with 
chiral extrapolations and with unquenching ap¬ 
pears to be much more favorable. Indeed, for 
Nf = 2, the coefficient of the chiral logarithm 
for Bb is in the range —0.2 < i(l — 3g^) < 0.4 
instead of 0.8 < |(1 -I- 3g^) < 1.9 as it is for /b, 
when 0.2 < g < 0.7. For g = 0.6, these coeffi¬ 
cients are 0.0 and 1.6, respectively. This milder 
behavior is confirmed by the unquenched calcu¬ 
lation of JLQCD [68,69], whose results for the 
light-quark-mass dependence of Bb and Bb^ are 
shown in Fig. 8. This mild light-quark-mass 
dependence is also observed in all quenched cal¬ 
culations [47, 48, 71-73] of which an example is 
shown in Fig. 8. Also evident from Fig. 8 is the 
fact the there is very little variation in going from 
Nf = 0 to Nf = 2, indicating that quenching ef- 

'•^The fact that the chiral behavior of the decay con¬ 
stants may be affected by the variation of lattice spac¬ 
ing and volume mentioned above does not imply that the 
B-parameters are (see discussion around Eq. (12)). 
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fects are small for Bg and Bg^ ■ 

The situation is less clear when it comes to the 
heavy-quark-mass behavior of Bg^. In [47,48], 
B parameters computed in the relativistic ap¬ 
proach for heavy-quark masses around that of 
the charm are extrapolated up to the 6, using 
HQET as a guide. Leading logs are accounted for 
in [48]. In [70], NLL corrections in the matching 
of QCD to HQET are included, and a calcula¬ 
tion of Bg in the static limit is used in conjunc¬ 
tion with the relativistic results of [47] to reach 
the 6-quark by interpolation. These three calcu¬ 
lations indicate that Bg increases with increas- 
ing Mg^. On the other hand, the NRQCD re¬ 
sults of [73], performed at quark masses above 
and around rrih, suggest an opposite behavior. 
This behavior is mild, however, and the different 
formulations agree at the physical 6-quark mass. 
Moreover, once systematics errors are included, 
the significance of these inconsistencies becomes 
marginal. 

The results of calculations of the NLO-RGI pa¬ 
rameter, Bg^°, are collected in Fig. 9, together 
with those for BgjBg. All methods give fully 
consistent results suggesting that residual system¬ 
atic errors cancel in the ratio of matrix elements, 
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Figure 9. Lattice results for Bff° and Bg^/Bg. 
Calculations are grouped by the approach used for 
the heavy quark. The results shown are: APE 00 
f47], LL 00 [48], APE 01 [70], BBS 98 [71], 
HASH 00 [72], JLQCD 01 [68], JLQCD 02 [73[. 
Also shown are my averages of the quenched re¬ 
sults for these quantities. 


of the conference. I use quenched averages as a 
starting point and the ratio of Nf = 2 to Nf = 0 
results to extrapolate these results to iV/ = 3, 
including an additional systematic error equal to 
the shift from Nf = 2 to Nf = 3 in the final 
result. I obtain 


^ ^ 3 {Bq\{bq)v_A{bq)v-A\Bq) 

8 (i3,lg7i"756l0)(0l67^75glR5) ’ 

used to define the R-parameters. This appears 
also to be true of quenching effects. 

The mild dependence on light valence-quark 
mass further guarantees that errors on BgjBg 
are small, of order 3%. The situation regard¬ 
ing the heavy-quark-mass dependence, however, 
should be clarified. This could be achieved 
by combining relativistic and HQET results ob¬ 
tained with non-perturbative renormalization in 
the continuum limit (see Sec. 2.3). 

4.3. Bq — mixing: summary 

I provide here summary numbers for the non- 
perturbative quantities relevant for the study of 
Bq — Hq-mixing. In performing averages, I have 
only kept results dated after 1998 in the quenched 
case {Nf = 0) and after 1999 for the unquenched 
case {Nf = 2) and omit results which had not 
made it into proceedings or papers at the time 


ffg = 203(27)t^oMeV, fg, = 238(31) Me^ 
Bff^° = 1.34(12), Bff^° = 1.34(12), 
/b\/^ = 235(33) j:°4MeV, 

^ Ib. = 276(38) MeV, J 


^ = 1.18(4)+^^ 
fs 


0 > 


Bg^ 

Bg 


= 1.00(3), 


C = = 1.18(4)1^2 




0 ’ 


where the second, asymmetric error is the one due 
to the uncertainty in the chiral extrapolation of 
fg discussed above. This error should be viewed 
as rough estimate of the range of possible devi¬ 
ations that can be induced by including chiral 
logs in this extrapolation. Thus, in UT analyses, 
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the value of ^ given above should be understood 
as ^ = 1.24(4)(6) and similarly for the other 
quantities affected by this error. The results were 
not presented in this way because the inclusion 
of chiral logs is still exploratory. Note that this 

discussion does not concern and /b, -8^^° 

whose residual chiral-log uncertainties appear to 
be small compared to other errors. 

Before closing, let me mention that two groups 
have calculated ^ directly from the ratio of AB = 
2 matrix elements [48,71]. The results obtained 
come out higher than the 1.18 given above, but 
are compatible within errors. 

4.4. K^-K^ —mixing and Bk 

The standard model, AS = 2 matrix element 
relevant for CP violation in — X°-mixing has 
been studied extensively on the lattice. A fea¬ 
ture common to all but one of the calculations 
reviewed below is the use of the quenched approx¬ 
imation. Furthermore, in these calculations, 
kaons are composed of mass-degenerate quarks 
and the final result is obtained by interpolation 
to a kaon of mass Mk, made up of two quarks 
of mass ~ nisi2. The procedure has the advan¬ 
tage that a chiral extrapolation to the mass of 
the down quark is avoided. The drawback is that 
S'[/(3)-breaking corrections ~ {ms — vridY Etre not 
included and must be accounted for in the sys¬ 
tematic error. 

4.4.1. New calculations of Bk 

This year three new quenched calculations of 
Bk have been performed. The first, with Wil¬ 
son type fermions, was presented at this confer¬ 
ence [3,75]. It is performed with high statistics 
at three different values of the lattice spacing in 
the range from 0.1 to 0.06 fermi, allowing for a 
continuum extrapolation. It implements a new 
method based on chiral Ward identities to elim¬ 
inate the spurious mixing with wrong chirality 
four-quark operators [76] as well as the more tra¬ 
ditional method of chiral subtractions. The au¬ 
thors find good agreement between their two pre¬ 
liminary results in the continuum limit. They 

'^^This result is entirely consistent with the recent phe¬ 
nomenological analysis of [74]. 

'^^The exception is the preliminary Nj = 3 result of [91]. 


renormalize the i3-parameter non-perturbatively 
in the RI/MOM scheme. 

The two other studies [77, 78] were performed 
using overlap quarks [79], a recent discretization 
of QCD which satisfies the Ginsparg-Wilson rela¬ 
tion [80] and thus exhibits an exact chiral-flavor 
symmetry, analogous to the one of continuum 
QCD, at finite lattice spacing [81]. This sym¬ 
metry implies that the AS' = 2 operator is sim¬ 
ple to construct on the lattice (unlike with stag¬ 
gered fermions) and does not mix with chirally 
dominant operators [82] (as it does with Wil¬ 
son fermions). Overlap fermions further guar¬ 
antee full O (a)-improvement. These calcula¬ 
tions are the first weak matrix element studies 
to use overlap quarks. They help validate this 
discretization as a useful phenomenological tool 
despite its present numerical cost. Such calcula¬ 
tions will eventually verify that the explicit break¬ 
ing of flavor (staggered fermions) or chiral (Wil¬ 
son fermions) symmetry is controlled in standard 
calculations. They will also provide a check of 
domain wall fermion calculations. Domain wall 
quarks [83,84] are a five-dimensional implemen¬ 
tation of Ginsparg-Wilson fermions. The chiral- 
flavor symmetry which they exhibit becomes ex¬ 
act in the limit of infinite fifth dimension. 

The first of these two overlap fermion calcula¬ 
tions [77] is performed in the quenched approx¬ 
imation at a single value of the lattice spacing 
{a~^ ~ 2 GeV) for five values of the light-quark 
mass in the range from TOs/ 2 to l.Sms- The B- 
parameter is renormalized non-perturbatively in 
the RI/MOM scheme. 

The second overlap calculation [78] is also per¬ 
formed in the quenched approximation, at two 
values of the lattice spacing (a“^ ~ 1.7 GeV and 
2.2 GeV), with a slightly different action and with 
five light-quark masses, ranging from 0.8ms to 
2.5ms. In the mass interval common to both cal¬ 
culations, agreement is excellent. However, the 
large quark masses considered in [78] imply that 
the physical point is reached through a somewhat 
dangerous (linear) extrapolation. The resulting 
R-parameter is renormalized at one-loop. Results 
at the two lattice spacings are in good agreement. 

Before summarizing the situation regarding 
Bk, I wish to mention that ALPHA [85] are per- 
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Figure 10. Lattice results for {2 G^) vs lat¬ 
tice spacing. The results are coded according to 
the quark action used: Ginsparg-Wilson (filled 
symbols), staggered or Kogut-Susskind (unfilled 
symbols) and Wilson (stick symbols). The contin¬ 
uum limit corresponds to a = 0. All but the refer¬ 
ence result of JLQCD [89] in that limit are shifted 
to the left for clarity. The results shown are: CP- 
PACS [87], RBC [88], GGHLR [77], MILC [78], 
JLQGD (KS) [89], KGS [90], OSU and iV/ = 3 
[91], KSGP [92], SPQR [75], JLQCD (Wilson) 
[93], UKQCD [94], BGS [95]. Not shown are the 
pioneering domain wall results of [96]. 

forming a high-statistics calculation of this B- 
parameter using non-perturbative renormaliza¬ 
tion and a modified version of Wilson fermions 
which goes under the name of twisted-mass QCD 
[86]. This approach circumvents the problem of 
mixing with wrong chirality operators in much 
the same way as the Ward identity method of [76], 
which it inspired. 

4.4.2. Summary of Bk determinations 

Results for R^°'^(2GeV) obtained by various 
groups using different approaches are plotted as 
a function of lattice spacing in Fig. 10. While re¬ 
sults differ at finite lattice spacing, those which 
are extrapolated to the continuum limit agree in 
this limit within one-and-some standard devia¬ 
tions. The overlap results of [77,78], obtained at 
finite lattice spacing, are also in good agreement 
with these continuum results, albeit with large 
statistical errors. Domain wall results are a bit 
on the low side. This difference should be un¬ 
derstood through further scrutiny of systematic 
errors in those calculations. 


The reference result is still the one obtained 
in the 1997 quenched staggered calculation of 
JLQCD [89], which is the most complete study of 
this quantity to date, despite its use of one-loop 
perturbation theory in matching the lattice oper¬ 
ators to their continuum counterparts. To their 
result one must add uncertainties due to quench¬ 
ing and the fact that the down and strange quarks 
are degenerate in the calculation. Sharpe [97], on 
the basis of quenched yPT and the preliminary 
“fV/ = 3” OSU results [91], suggests an enhance¬ 
ment factor of 1.05±0.15 to “unquench” quenched 
results for Bk ■ He advocates an additional factor 
1.05 ± 0.05 to compensate for the effect of work¬ 
ing with degenerate quark masses. Here, I choose 
to keep his estimate of errors, but not to use his 
enhancement factors so as not to modify the cen¬ 
tral value on the basis of information which is 
provided by yPT estimates and preliminary lat¬ 
tice results. Thus, I quote 

A \ 

H^°'^(2GeV) = 0.628(42)(99) 

^ ^ Bf^° = 0.86(6) (14) , _^ 

where B^° is the two-loop RGI R-parameter 
obtained from GeV) with Nf = 3 and 

as(2GeV) = 0.3. 

This is the same result as was presented at Lat¬ 
tice 2000 [7]. What is really needed now are high 
statistics unquenched calculations to reduce the 
dominant 15% quenching error. Without such an 
improvement in accuracy, indirect CP violation 
in the neutral kaon system will no longer provide 
a significant constraint on UT fits [44]. 

4.5. Semileptonic decays 
4.5.1. \Vcb\ from B D*£v 

The CKM element \ Vcb\ plays an important role 
in constraining the UT, for it determines the tri¬ 
angle’s base but also appears raised to the fourth 
power in the constraint on the summit provided 
by Cif ■ jbcfej must thus be determined as precisely 
as possible. 

One way to measure \Vcb\ accurately is by ex¬ 
trapolating the differential decay rate for B 
D*£v decays. 
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to the zero-recoil point w = vb ■ f_D* = 1, where 
vb and vb* are the four-velocities of the ini¬ 
tial and final mesons [98]. The benefit of such 
an extrapolation is that heavy quark symmetry 
[99, 100] requires that the form factor To* {w) 
must be normalized to 1 at ic = 1, up to small ra¬ 
diative and small order ^Iml^ corrections [101]. 
Nevertheless, a few percent measurement of |14h| 
requires a reliable determination of (1) — 1. 
While the radiative corrections have been calcu¬ 
lated to two-loops in perturbation theory [102], 
the long-distance power corrections require a non- 
perturbative treatment. 

Through a clever use of double ratios of matrix 
elements for D^*\ B^*'> D^*'>, B^*'>-like transi¬ 

tions, the authors of [103] have been able to ob¬ 
tain a statistically significant signal for the lead¬ 
ing non-perturbative corrections of order 1 /ml ^ 
and some of the sub-leading ones proportional to 
three factors of the inverse heavy-quark masses. 
The result of their quenched calculation, per¬ 
formed at three values of the lattice spacing, was 
presented at this conference [4,103]: 

.Fzr.(l) = 0.913l?^t" . 

This important result, which requires excellent 
control of statistical and systematic errors, should 
be verified by other lattice groups. 

4.6. B irii/ and \Vub\ 

In heavy-to-light semileptonic decays, the light, 
final state hadron can have momenta as large as 
|p[ ~ toq/ 2 in the parent rest frame, where mg 
is the mass of the heavy quark. For mg = mb, 
and on present day lattices, such momenta would 
lead to uncontrollably large discretization effects, 
proportional to powers of a\p\. Therefore, at 
present, only a limited kinematical range below 
the zero-recoil point can be reached without ex¬ 
trapolation. Even so, lattice calculations are use¬ 
ful, for the relevant matrix elements are not nor¬ 
malized by heavy quark symmetry as they are in 
B decays. Furthermore, experiment is 

beginning to measure the corresponding differen¬ 
tial rates in the kinematical region accessible to 
lattice calculations (see e.g. [104]), thus allowing 
model-independent determinations of |I4b|- 
The matrix element relevant for B^ —>■ 



form factors F+{q^) and Fq^q^) (figure from [3]). 
These results are APE [105], UKQCD [106], 
FNAL [107], JLQCD [108]. Also shown is a fit 
to the parametrization of [109] and the light-cone 
sumrule results of [110] (see also [111]). 

decays is 

{lT-{p')\V^^\B\p)) = 

+ {p + p' , 

where q = p — p' and = Iry^u. There are four 
quenched calculations of this matrix element, us¬ 
ing either relativistic, Fermilab or NRQCD heavy 
quarks. The form factors Fj^{qY and Fo(q^) ob¬ 
tained by the different collaborations are plotted 
as functions of q^ in Fig. 11. Excellent agreement 
is found for Fl|_(g^) which determines the rate in 
the limit of vanishing lepton mass. This agree¬ 
ment is certainly not trivial given the use of dif¬ 
ferent heavy-quark approaches and the variations 
in the rather intricate analyses leading to the de¬ 
termination of these form factors. The error on 
this form factor is typically of order 20%. 

The lattice results cover roughly the upper half 
of the kinematically accessible range. The easi¬ 
est way to extrapolate these results to lower val¬ 
ues of q^ is to use a parametrization such as the 
one presented in [109], which incorporates most of 
the known constraints on the form factors. Such 
an extrapolation is shown in Fig. 11 and gives 
excellent agreement with light-cone sumrule re¬ 
sults [110,111]. This extrapolation introduces. 
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however, a model dependence. To eliminate this 
model dependence, one can make use of the dis¬ 
persive bound techniques proposed in [112,113]. 
In fact, the results of the analysis of [113] could 
certainly be improved with the lattice results pre¬ 
sented here. 

Given the good agreement amongst the dif¬ 
ferent quenched calculations of F+{q^), what is 
needed to reduce the error on the extraction of 
\Vub\ are unquenched calculations and more work 
on methods for enlarging the kinematical reach of 
lattice calculations. 

5. Outlook 

A number of quantities of central importance to 
particle physics are currently being studied with 
lattice QCD, many of which could unfortunately 
not be presented here. For those that were, em¬ 
phasis is, in most cases, on the reduction of sys¬ 
tematic errors, such as quenching, discretization 
errors, etc. 

Because they have not yet led to phenomeno¬ 
logical breakthroughs, I did not discuss in any 
detail the major advances of the last few years 
associated with the formulation and implementa¬ 
tion of fermionic discretizations which preserve 
an exact continuum-like chiral-flavor symmetry 
at finite lattice spacing. These discretizations are 
known as Ginsparg-Wilson fermions [80], of which 
overlap [79], domain wall [83] and fixed-point 
fermions [114] are explicit realizations. They have 
opened new possibilities for the calculation of 
weak matrix elements, and in particular those as¬ 
sociated with the A/ = 1/2 rule and direct CP 
violation in K tttt decays. On the analytical 
side, the renormalization of the weak operators 
involved has been studied with domain wall [115] 
and overlap fermions [116,117]. On the numerical 
side, a first round of investigations has been per¬ 
formed with domain wall fermions [88,118], laying 
the foundations for future studies. 

The exceptional chiral properties of these new 
fermion formulations has also allowed the investi¬ 
gation of an unexplored regime of QCD in which 
the correlation length of the pion field is much 
larger than the linear extent of spacetime, the e- 
regime of [119]. This investigation has lead to 


exploratory calculations of one of the low-energy 
constants (LEG) of the strong chiral lagrangian 
in the quenched approximation [120-123]. It is 
possible, in principle, to generalize this approach 
to extract the LECs of the weak chiral lagrangian 
and a numerical investigation is under way [124]. 

More generally, the range of approaches and of 
quantities studied in lattice QCD is constantly ex¬ 
panding and there should be many more exciting 
results to present at ICHEP 2004. 
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